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The long-distance contributions to K — > 2tt amplitudes can be pinned 
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strong S-wave rescattering of the two final pions generates sizeable chiral 
loop corrections, which have an important impact on the direct CP vio- 
lation ratio e'/s Including all large logarithmic corrections, both at 
short and long distances, the Standard Model Prediction for this observ- 
able is found to be Re(e'/e) = (1.7 ±0.9) • 10~ 3 , in good agreement 
with the most recent experimental measurements. A better estimate of 
the strange quark mass could reduce the theoretical uncertainty to 30%. 
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1 Introduction 



The CP-violating ratio e' / e constitutes a fundamental test for our understanding 
of flavour-changing phenomena within the Standard Model framework. The experi- 
mental status has been clarified by the recent KTEV |§, Re (e'/e) = (28.0±4.1)-10~ 4 , 
and NA48 ||], Re {e'/e) = (15.3 ± 2.6) • 10~ 4 , measurements, which provide clear ev- 
idence for a non-zero value and, therefore, the existence of direct CP violation. The 
present world average is H[|,|5],|6|] , 

Re (e'/e) = (18.0 ±2.0) • 1(T 4 , (x7 n df = 10.8/3). (1) 

The theoretical prediction has been rather controversial since different groups, us- 
ing different models or approximations, have obtained different results |7|,§,[j[ [L0[[Ll , 12] . 



Although there was no universal agreement on the e'/e value predicted by the Stan- 
dard Model, it has been often claimed that it is too small, failing to reproduce the 
experimental world average by at least a factor of two. This claim has generated a 
very intense theoretical activity, searching for new sources of CP violation beyond 
the Standard Model pj. 

It has been pointed out recently that the theoretical short-distance evalua- 
tions of e'/e had overlooked the important role of final-state interactions (FSI) in 
K — > mi decays. Although it has been known for more than a decade that the 
rescattering of the two final pions induces a large correction to the isospin-zero decay 
amplitude, this effect was not taken properly into account in the theoretical predic- 
tions. From the measured ti-ti phase shifts one can easily infer |IJ that FSI generate a 
strong enhancement of the e'/e prediction, by roughly the needed factor of two. This 
large correction is associated with infrared chiral logarithms involving the pion mass, 
which can be rigorously analyzed with standard Chiral Perturbation Theory (xPT) 
techniques |]T^,|T5],|l6[ . A very detailed analysis, including all large logarithmic correc- 



tions both at short and long distances, has been presented in ref. The resulting 
Standard Model prediction [0, 

Re (e'/e) = (17 ±9) ■ 10~ 4 , (2) 

is in good agreement with the most recent measurements. 

The following sections present a brief overview of the most important ingredients 
entering the theoretical prediction of e'/ e : 

1. A short-distance calculation at the electroweak scale and its renormalization- 
group evolution to the three-flavour theory (/i < m c ), which sums the leading 
and next-to-leading ultraviolet logarithms. 

2. The matching to the ;^PT description, which so far has been done at leading 
order in the 1/Nc expansion. 

3. Chiral loop corrections, which induce large infrared logarithms related to FSI. 
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2 Theoretical framework 



In terms of the K — > tttt isospin amplitudes, Ai = Aj e Sl (I = 0, 2) 

~lm(A 2 ) Im(A ) 



e' to 



e 1 * 



V2\e 



Re(A 2 ) Re(A ) 



(3) 



Owing to the well-known "A/ = 1/2 rule", e'/e is suppressed by the ratio uj = 
Re(A 2 ) /Re(A ) « 1/22. The strong S-wave rescattering of the two final pions gen- 
erates a large phase-shift difference between the two isospin amplitudes, making the 
phases of e' and e nearly equal. Thus, 

$«5 2 -5 + ^~0. (4) 

The CP-conserving amplitudes Re(Aj), their ratio u and e are usually set to their 
experimentally determined values. A theoretical calculation is then only needed for 
the quantities Im(Aj). 

One starts above the electroweak scale where the flavour-changing process, in 
terms of quarks, leptons and gauge bosons, can be analyzed within the usual gauge- 
coupling perturbative expansion in a rather straightforward way. Since Mz is much 
larger than the long-distance hadronic scale Mk, there are large short-distance loga- 
rithmic contributions which can be summed up using the Operator Product Expansion 
(OPE) |l7j and the renormalization group. The proper way to proceed makes use of 
modern Effective Field Theory (EFT) techniques |18| . 



The renormalization group is used to evolve down in energy from the electroweak 
scale, where the top quark and the Z and W ± bosons are integrated out. That means 
that one changes to a different EFT where those heavy particles are no longer explicit 
degrees of freedom. The new Lagrangian contains a tower of operators constructed 
with the light fields only, which scale as powers of 1/Mz- The information on the 
heavy fields is hidden in their (Wilson) coefficients, which are fixed by "matching" 
the high- and low-energy theories at the point fi = Mz- One follows the evolution 
further to lower energies, using the EFT renormalization group equations, until a new 
particle threshold is encountered. Then, the whole procedure of integrating the new 
heavy scale and matching to another EFT starts again. 

One proceeds down to scales \x < m c and gets finally an effective AS" = 1 La- 
grangian, defined in the three-flavour theory Hl9| , |20| . 



fi 10 

#/f 1 = — 7=Kd C E , (5) 

V2 i=l 

which is a sum of local four-fermion operators Qi, constructed with the light degrees 
of freedom, modulated by Wilson coefficients Cj(/x) which are functions of the heavy 
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Figure 1: Evolution from Mg to Mj(. 



masses. We have explicitly factored out the Fermi coupling G7p and the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements containing the usual Cabibbo sup- 
pression of K decays. The unitarity of the CKM matrix allows to write 

Ci{p) = Zi{p) +r yi{n) , (6) 

where r = —V td V t */V ud V* s . Only the ^ components are needed to determine the 
CP-violating decay amplitudes. The overall renormalization scale /i separates the 
short- (M > fi) and long- (m < fx) distance contributions, which are contained in 
Ci(n) and Qi, respectively. The physical amplitudes are of course independent of \i. 

Our knowledge of AS* = 1 transitions has improved qualitatively in recent years, 
thanks to the completion of the next-to-leading logarithmic order calculation of the 
Wilson coefficients pT| , p2| . All gluonic corrections of 0(a™t n ) and 0(a™ +1 t n ) are 
known, where t = In (Mi/M 2 ) refers to the logarithm of any ratio of heavy mass 
scales M 1 ,M 2 > fi. Moreover, the full m t /M w dependence (at lowest order in a s ) is 
taken into account. 

In order to predict physical amplitudes, however, one is still confronted with the 
calculation of hadronic matrix elements of the four-quark operators. This is a very 
difficult problem, which so far remains unsolved. Those matrix elements are usually 
parameterized in terms of the so-called bag parameters Bi, which measure them in 
units of their vacuum insertion approximation values. 
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To a very good approximation, the Standard Model prediction for e'/e can be 
written (up to global factors) as |7|] 



n 



IB 



0.4i4 3/2) 



(7) 



Thus, only two operators are numerically relevant: the QCD penguin operator Q 6 
(AI = 1/2) governs lm(Ao), while lm(A2) (AI = 3/2) is dominated by the elec- 
troweak penguin operator Qg. The parameter 



n 



IB 



1 Im(A 2 ) IB 
uj Im(v4o) 



takes into account isospin-breaking corrections, which get enhanced by the large 
factor \/u. 

The isospin-breaking correction coming from TT°-r] mixing was originally estimated 
to be = 0.25 p3| , p4| . Together with the usual ansatz Bi ~ 1, this produces a 
large numerical cancellation in eq. (|7[) leading to low values of e'/e around 7- 10~ 4 . A 
recent improved calculation of ir -r] mixing at 0(p A ) in has found the result 



IB 



0.16 ±0.03. 



(9) 



This smaller number, slightly increases the naive estimate of e'/e. 



3 Chiral Perturbation Theory 

Below the resonance region one can use global symmetry considerations to define 
another EFT in terms of the QCD Goldstone bosons (7r, K, if). The xPT formulation 
of the Standard Model [|n],|l^,|T(J describes the pseudoscalar-octet dynamics, through 
a perturbative expansion in powers of momenta and quark masses over the chiral 
symmetry breaking scale (A x ~ 1 GeV). 

Chiral symmetry fixes the allowed %PT operators. At lowest order in the chiral 
expansion, the most general effective bosonic Lagrangian with the same SU(3)l ® 
SU (3)r transformation properties as the short-distance Lagrangian (||) contains three 
terms, transforming as (8^,1^), (27^,1^) and (8^,8^), respectively. Their corre- 
sponding chiral couplings are denoted by g$, g^j and g ew . 

The tree-level K — > tttt amplitudes generated by the lowest-order %PT Lagrangian 
do not contain any strong phases: 

Ao = -% V ud V: s V2U {(gs + l g 27 ) (M* - M 2 % ) - 2 - ft e 2 g 8 g ew } , 

A ^ = -% V udKs I U {s 927 {M 2 K — Ml) — 3 fl e 2 g 8 g ew } . (10) 
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Taking the measured phase shifts into account, the moduli of gs and #27 can be 
extracted from the CP-conserving K — > 2tt decay rates; one gets |26j \g 8 \ ^5.1 
and I #27 1 ~ 0.29. The huge difference between these two couplings shows the well- 
known enhancement of octet |A/| = 1/2 transitions. The g ew term is the low-energy 
realization of the electroweak penguin operator. 

The isospin amplitudes Ai have been computed up to next-to-leading order in 
the chiral expansion P, p7| , p^j29|j30|j3T[| . The only remaining problem is the calculation 
of the ;\TT couplings from the effective short-distance Lagrangian (||), which requires 
to perform the matching between the two EFTs. This can be easily done in the large- 
Nc limit of QCD [fffil , because in this limit the four-quark operators factorize into 
currents which have well-known chiral realizations. The local 0(p A ) contributions to 
the amplitudes Ai can be easily included in eqs. (|10"D , through effective correction 
factors AqA^ to the lowest-order g^ contributions. At leading order in 1/Nq, one 
gets O: 



I c x <ji) + 1 cm + cm - 16 l 5 cm 
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where f K ' K (Ml) « 1 +4L 5 Ml/ fl is the Ktt scalar form factor at the pion mass scale, 



L 5 is a coupling of the 0{p ) strong chiral Lagrangian and m q 



m,, 



rrid- In the 



limit N c -> 00, Lf = \fl - l) / {M 2 K - Ml) « 2.1 • 1(T 3 and ~ 1-02. 

These results are equivalent to the usual large-iVc evaluations of the Bi factors. 
In particular, for e'/e where only the imaginary part of the gn couplings matter [i.e. 



Im(Cj)] they amount to B 



(3/2) 



B, 



(1/2) 



1. Therefore, up to minor variations on 



some input parameters, the corresponding e'/e prediction reproduces the published 
results of the Munich M and Rome M groups. 
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The large-iVc limit is only applied to the matching between the 3-flavour quark 
theory and xPT, as indicated in Figure [I]. The evolution from the electroweak scale 
down to fi < m c has to be done without any unnecessary expansion in powers of 
1/Nc', otherwise, one would miss large corrections of the form In (M/m), with 
M ^> m two widely separated scales ]33]]. Thus, the Wilson coefficients contain the 
full \i dependence. 

The large- N c factorization of the four-quark operators Q t (i ^ 6, 8) does not 
provide any scale dependence. Since the anomalous dimensions of these operators 
vanish when Nc — > oo []33|1 , a very important ingredient is lost in this limit |34 . 
To achieve a reliable expansion in powers of 1/Nc, one needs to go to the next 
order where this physics is captured p4|,35[ . This is the reason why the study of 
the AI = 1/2 rule has proven to be so difficult. Fortunately, these operators are 
numerically suppressed in the e'/s prediction and their contributions can be easily 
corrected with the information provided by the measured CP-conserving rates pjjTj. 

The only anomalous dimensions which survive when Nc — > oo are precisely the 
ones corresponding to Qq and Q$ p4| , p3| . One can then expect that the matrix 
elements of these two operators are well approximated by this limit [^,^,^| . These 
operators factorize into colour-singlet scalar and pseudoscalar currents, which are 
/i dependent. This generates the factors (qq)(fi) ~ — f^M\j{m s + m q )(fi), which 
exactly cancel the \i dependence of C%${p) at large Nc P^ , |33| , p^ , |35| , p6| , p7| . It remains 
of course a dependence at next-to-leading order. 

Therefore, while there are large 1/Nc corrections to Re(<7/) [P5[j , the large-iVc 
limit is expected to give a good estimate of Im(<7j). 



4 Chiral loop corrections 

The lowest-order calculation does not provide any strong phases 8j. Those phases 
originate in the final rescattering of the two pions and, therefore, are generated by 
chiral loops which are of higher order in both the momentum and 1/Nc expansions. 
Analyticity and unitarity require the presence of a corresponding dispersive FSI effect 
in the moduli of the isospin amplitudes. Since the S-wave strong phases are quite 
large, specially in the isospin-zero case, one should expect large higher-order unitarity 
corrections. 

The size of the FSI effect can be calculated at one loop in xPT. The dominant one- 
loop correction to the octet amplitude comes indeed from the elastic soft rescattering 
of the two pions in the final state. The existing one-loop analyses P J27|J28|] show that 
pion loop diagrams provide an important enhancement of the Aq amplitude, implying 
a corresponding reduction of the phenomenologically fitted value of \g%\. This chiral 
loop correction destroys the accidental numerical cancellation in eq. (0), generating 
a sizeable enhancement of the e'/e prediction M. 
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Let us decompose the isospin amplitudes in their different chiral components as 

— M.Q 

the form 



Aq = A^ + Aq + A^ w ^ and A2 = A^ + A^^ ■ Moreover, we can write them in 



A (R) = A (R)oo x C (R) ^ (12) 

where Aj°° are the large- Nc results obtained in the previous section. The correction 
factors C\ = 1 + A L A} contain the chiral loop contributions A L A) that we are 
interested in. Their complete analytical expressions at one loop in %PT have been 
given in ref. J|, where the following numerical values have been obtained: 

Cj 8) = 1.27 ±0.05 + 0.46 i, 

Cj 27) = 2.0 ±0.7 + 0.46 i, 

C^ w) = 1.27 + 0.05 + 0.46 i, 

Cf 7) = 0.96 + 0.05 - 0.20 i, 

C { 2 ew) = 0.50 + 0.24- 0.20 i. (13) 

The central values have been evaluated at the chiral renormalization scale v = M p . 
To estimate the corresponding uncertainties we have allowed the scale v to change 
between 0.6 and 1 GeV. The scale dependence is only present in the dispersive con- 
tributions and should cancel with the corresponding v dependence of the local %PT 
counterterms. However, this dependence is next-to- leading in 1/Nc and, therefore, is 
not included in the large- Nc determination of the chiral couplings. The sensitivity of 
the results to the scale v gives a good estimate of those missing contributions. Notice 
that all amplitudes with a given isospin get the same absorptive contribution, as it 
should since they have identical strong phase shifts. 

The numerical corrections to the 27-plet amplitudes do not have much pheno- 
menological interest for CP-violating observables, because Im(g27) = 0. Remember 
that the CP-conserving amplitudes Ke(Ai) are set to their experimentally deter- 
mined values. What is relevant for the e'/e prediction is the 35% enhancement of the 
isoscalar octet amplitude Im^g 8 -*] and the 46% reduction of ImfAjT"'']. Just looking 
to the simplified formula (0), one realizes immediately the obvious impact of these 
one-loop chiral corrections. 

A complete 0(p 4 ) calculation [^,[^J of the isospin-breaking parameter Qjb is 
not yet available. The value 0.16 quoted in eq. only accounts for the contribution 
from 7r°-?7 mixing |25| and should be corrected by the effect of chiral loops. Since 
{C^l ~ 0.98 + 0.05, one does not expect any large correction of lm(A 2 )iB, while we 
know that Im[Ao 8 ^] gets enhanced by a factor 1.35. Taking this into account, one gets 
the corrected value 

/7(27) 

Zl 
M8) 
u 



■7T°ri 

IB ~ "/_B 



0.12 + 0.05, (14) 
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where the quoted error is an educated theoretical guess. This value agrees with the 
result Qjb = 0.08 ± 0.05 ± 0.01, obtained in ref. |38| by using three different models 
P j3l| , |35| , |39l , ^0| , ^l| to estimate the relevant 0(p 4 ) chiral couplings. 



The sensitivity to higher-order chiral loop corrections has been investigated in 
ref. p| through an Omnes exponentiation of the dominant pion loops 0, using the 
experimental tttt phase shifts. The standard one-loop results and the Omnes 
calculation agree within errors, indicating a good convergence of the chiral expansion. 



5 Numerical results and discussion 

The infrared effect of chiral loops generates an important enhancement of the 
isoscalar K — > tttt amplitude. This effect gets amplified in the prediction of s'/s, 
because at lowest order (in both 1/Nc and the chiral expansion) there is an accidental 
numerical cancellation between the 1 = and 1 = 2 contributions. Since the chiral 
loop corrections destroy this cancellation, the final result for e'/ e is dominated by the 
isoscalar amplitude. Thus, the Standard Model prediction for e' / e is finally governed 
by the matrix element of the gluonic penguin operator Q e . 

A detailed numerical analysis has been provided in ref. 0. The short-distance 
Wilson coefficients have been evaluated at the scale \l = 1 GeV. Their associated 
uncertainties have been estimated through the sensitivity to changes of \i in the 
range M p < // < m c and to the choice of 75 scheme. Since the most important 
a s corrections appear at the low-energy scale //, the strong coupling has been fixed 



at the r mass, where it is known |42|] with about a few percent level of accuracy: 
a s {m T ) = 0.345 ± 0.020. The values of a s at the other needed scales can be deduced 
through the standard renormalization group evolution. 

Taking the experimental value of e, the CP-violating ratio e'/e is proportional to 



the CKM factor Im(^*l^) = (1.2 ± 0.2) • 10~ 4 J43|. This number is sensitive to the 
input values of several non-perturbative hadronic parameters adopted in the usual 
unitarity triangle analysis; thus, it is subject to large theoretical uncertainties which 
are difficult to quantify [Q. Using instead the theoretical prediction of e, this CKM 
factor drops out from the ratio e'/e] the sensitivity to hadronic inputs is then reduced 
to the explicit remaining dependence on the AS* = 2 scale-invariant bag parameter 
Bk- In the large- Nq limit, B^ = 3/4. We have performed the two types of numerical 
analysis, obtaining consistent results. This allows us to estimate better the theoretical 
uncertainties, since the two analyses have different sensitivity to hadronic inputs. 
The final result quoted in ref. is: 

Re (e'/e) = (l.7 ± 0.2 ±gj ± 0.5) • 10~ 3 = (1.7 ± 0.9) • 10^ 3 . (15) 

The first error comes from the short-distance evaluation of Wilson coefficients and 
the choice of low-energy matching scale /1. The uncertainty coming from vary- 
ing the strange quark mass in the interval (m s + m g )(lGeV) = 156 ± 25MeV 



45| , ^6j^7| , ^8i ^9 lp0i^l| is indicated by the second error. The most critical step is the 
matching between the short- and long-distance descriptions. We have performed this 
matching at leading order in the 1/Nc expansion, where the result is known to 0(p 4 ) 
and 0(e 2 p 2 ) in x?T. This can be expected to provide a good approximation to the 
matrix elements of the leading Qq and Q& operators. Since all ultraviolet and infrared 
logarithms have been resummed, our educated guess for the theoretical uncertainty 
associated with l/N c corrections is ~ 30% (third error). 

Thus, a better determination of the strange quark mass would allow to reduce the 
uncertainty to the 30% level. In order to get a more accurate prediction, it would 
be necessary to have a good analysis of next-to-leading 1/Nc corrections. This is 
a very difficult task, but progress in this direction can be expected in the next few 
years 0,0,HJ52|H,|4 • 



To summarize, using a well defined computational scheme, it has been possible to 
pin down the value of e' / e with an acceptable accuracy. Within the present uncertain- 
ties, the resulting Standard Model theoretical prediction flT5|) is in good agreement 
with the measured experimental value (|l|), without any need to invocate a new physics 
source of CP violation. 
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